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Abstract 

We show that, under certain smoothness conditions, a Brownian martingale at a fixed time can be 
represented as an exponential of its value at a later time. The time-dependent generator of this exponential 
operator is equal to one half times the Malliavin derivative. This result can also be seen as a generalization 
of the semi-group theory of parabolic partial differential equations to the parabolic path-dependent partial 
differential equations introduced by Dupire (2009) and Cont and Founie (2011). The exponential operator 
. can be calculated explicitly in a series expansion, which resembles the Dyson series of quantum mechanics. 

0^ ' Our continuous-time martingale representation result is proved by a passage to the limit of a special case 

of a backward Taylor expansion of an approximating discrete-time martingale. The latter expansion can 
also be used for numerical calculations. 

Keywords: Continuous Martingales, Malliavin Calculus, Probabilistic Methods for Partial Differen- 
tial Equations, Path-Dependent Partial Differential Equations. 
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1 Introduction 

in 

■ The problem of representing Brownian martingales has a long and distinguished history. Dambis (1965) and 

Dubins-Schwarz (1965) showed that continuous martingales can be represented as time-changed Brownian 
1^ ■ motions. Doob (1953), Wiener and Ito developed what is often called Ito's martingale representation theorem: 

every local Brownian martingale has a version which can be written as an Ito integral plus a constant. In 
. this article, we consider martingales which are conditional expectations of a J-V-measurable random variable 

F. When the random variable F is Malliavin differentiable, the Clark-Ocone formula (Clark (1970), and 
Ocone (1984)) states that the integrand in Ito's martingale representation theorem is equal to the conditional 
expectation of the Malliavin derivative of F. We consider the less general problem of "infinitely smooth" 
martingales, namely martingales which are conditional expectations of a J^-measurable random variable F, 
which is infinitely differentiable in the sense of Malliavin. We show that a Brownian martingale at a fixed 
time can be represented as an exponential of its value at a later time. The time-dependent generator of 
this exponential operator is equal to one half times the Malliavin derivative. While smoothness is a severe 
limitation to our result, our representation formula opens the way to new numerical schemes, and potentially 
some analytical asymptotic calculations. The exponential operator can be calculated explicitly in a series 
expansion, which resembles the Dyson series of quantum mechanics. There are two main differences between 
our martingale representation and the Dyson formula for the initial value problem in quantum mechanics. 
First, in the case of martingales, time flows backward. Secondly, the time-evolution operator is equal to 
one half of the second-order Malliavin derivative, while for the initial value problem in quantum mechanics 
the time-evolution operator is equal to —2ni times the time-dependent Hamiltonian divided by the Planck 
constant. 

Our continuous-time martingale representation result is proved by a passage to the limit of a special 
case of the backward Taylor expansion (BTE) of an approximating discrete-time martingale. We introduced 
(without formal proof) the BTE in Schellhorn and Morris (2009), and applied it to price American options 
numerically. The idea in that paper was to use the BTE to approximate, over one time-step, the conditional 
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expectation of the option value at the next time-step. While not ideal to price American options because of 
the lack of differentiability of the payoff, the BTE is better suited to the numerical calculation of the solution 
of smooth backward stochastic differential equations (BSDE). In a related paper, Hu, Nualart, and Song 
(2011) introduce a numerical scheme to solve a BSDE with drift using Malliavin calculus. Their scheme can 
be viewed as a Taylor expansion carried out until the first order. Our BTE can be seen as a generalization to 
higher order of that idea, where the Malliavin derivative(s) is (are) calculated at the future time-step rather 
than at the current time-step. 

In the second part of the paper, we consider the path-dependent partial differential equation (PPDE) 
introduced by Dupire (2009). By the functional Feynman-Kac formula introduced by Dupire(2009) and Peng 
and Wang (2011), the solution of the PPDE is directly related to the martingale solving the corresponding 
BSDE. We can thus characterize the solution of the PPDE as an exponential of the terminal condition. This 
result belongs to the classical domain of integration of the evolution equation (see e.g. Yosida (1978), chapter 
14). In the smooth case, it generalizes in a natural way the semi-group theory of integration of the classical 
diffusion equation (see e.g., McOwcn (1996), chapter 9.2). 

The structure of this paper is the following. We first expose the discrete-time result, namely the Backward 
Taylor Expansion (BTE) for discrete functionals, and then prove our main result, namely our martingale 
representation theorem. Two explicit examples are given, which show the usefulness of the Dyson series in 
analytic calculations. In the second part of the paper we review the definition of the PPDE and provide an 
explicit exponential formula for its solution. The PPDE notation, which is not necessary to understand the 
first part, is introduced in the second part. The proofs of the theorems, as well as the notation necessary for 
the proofs, are relegated to the appendix. 



2 Martingale Representation 

The uncertainty is described by the filtered probability space (O, F, {F t }, P). Here {F t } is the filtration 
generated by Brownian motion W on K. Most results can be easily generalized to Brownian motion on R d . 
We denote the Malliavin derivative of order I of F at time t by D\F . We call Doo([0,r]) the set of random 
variables which are infinitely Malliavin differentiable and J^-measurable. Given u> <G Q, we denote by w t (u>) 
be the path that "freezes" Brownian motion after time t: 

ti w f W(s,uj) .. s<t 
*WM) = { W \tJ * t <7<T 

As is conventional for regular derivatives, the notation D l t F (w* (uo)) refers to the value of the Malliavin 
derivative of F along scenario cj*(cj), and not to the value of the composition of Fow* evaluated at lo. For 
instance \D\ VF 2 (w*(u;)) = l[s < T]W(t,u). 



2.1 Backward Taylor Expansion 

THEOREM 1 Let F £ B>oo([0, T])be a a {W (A), ..,W(T)} -measurable random variable. Suppose that 
D| m+1 - )A F = for I > L. Let j(m, L) be given by 7(m, L) = 1 if L = and, otherwise by: 

7 (m,L) = l{Leven}^-j _ - g 7 ( ,Z) ^-^ (1) 

Then 

L 

E[F\F mA ] = ]T 7 (m, l)E[D\ m+1)A F\F {{m+1)A) } (2) 

1=0 

We now consider functionals F which Malliavin derivatives do not vanish. The main application of a Taylor 
series comes from truncating it. We now proceed to estimate the truncation error. Let F^ be an approxima- 
tion of E[F\F m A] obtained by supposing that F has order L: 
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L 

F m = ^2'r(m,l)E[Dl m+1)A F\J : ^ m+1)A) } 

1=0 

THEOREM 2 Let F e D oo ([0,T])&e a a{W(A), .., W(T)} -measurable random variable. The mean square 
truncation error is: 

E[(E[F\F rnA ] - F r L n ) 2 \F mA ] = 0(A L+1 ) (3) 
Combining theorems 1 and 2, we arrive at the backward Taylor expansion. 

THEOREM 3 Let F E D^QO, T])be a a {W( A), ..,W(T)} -measurable random variable. Then 

OO 

E[F\F mA ] = Z)£[£>( m+1)A F|.F ((m+1)A) ] (4) 

1=0 

Applying (j4| recursively, one obtains the following corollary. 
COROLLARY 

oo oo M 

E[F\T mA \ = J2 "E II ^^J^^da-^a^ (5) 

3m+l=0 Ju=o i=m+l 

Observations: A non-intuitive feature of the Backward Taylor expansion is that any path can be chosen 
to approximate conditional expectations backward. Suppose that the paths of Brownian motion are fixed 
in advance, like in regression-based algorithms to calculate American options. The BTE can be used as a 
type of "control variate" to speed up the convergence of the backward induction, as we shall make precise in 
another article. If the paths are not chosen in advance, an interesting equation similar to ([7]) emerges when 
we choose the paths: 

W((m + 1)A) = W(mA) + iVE 
In the next subsection we will choose the "certainty-equivalent" paths: 

W((m+1)A) = W(mA) 

to derive our main result. For numerical applications, it is crucial to choose a low order of expansion L, in 
order to keep the number of calculations in (|5|) from growing too fast. One could then imagine a scheme 
where, at each step, the "optimal path" is chosen so as to minimize the global truncation error in We 
leave all these considerations for future research. 

2.2 Exponential Formula 

For esthetical reasons we introduce a "chronological operator". In this we follow Zeidler (2006). Let (H(t)) 
be a collection of operators. The chronological operator T is defined by 

r(H(t l )H{t 2 )..H{t n )) := H(t v )H{t 2 >)...H{t n ,) 

where ty, ..,t n > is a permutation of t\, ..,t n such that ty >t 2 ' > •• > t n i. 
Example 

It is showed in Zeidler (2006) p. 44-45 that 

/ / 2 fr(*i)£r(* 2 )diidta = \ f [ TiHit^H^dhdh 
Jo Jo ^ Jo Jo 

This will be the only property of the chronological operator we will use in this article. 
Definition: the exponential operator of a time-dependent generator H is: 
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rexp(/ H(s)ds) = J2 ■ T(H(r 1 ),..H(T k ))dT 1 ..dT k (6) 
Jt k=0 Jt Jt 

In quantum field theory, the series on the right handsidc of ([6J is called a Dyson series. 
THEOREM 3 Suppose F <= £>«,([(), T]). Lei Af(i) = B^l-Ft] /or i < T. Then, in L 2 {P) : 

M{t,w)=TeM\J Dlds)F{w\u)) (7) 

The importance of the exponential formula ([7]) stems from the Dyson series representation (j6j), which we 
rewrite hereafter in a more convenient way: 

M(t,Lj)=F(u t (uj)) + ±£ D 2 s F{u\u))ds+\j^ jf D 2 sl D 2 S2 F{uj t {uj))ds2ds 1 + ... 

It can be used for either numerical calculations (which we have not tried yet) of for analytical calculations, 
as we show in the next subsection. 

2.2.1 Solution of Some Problems by Dyson Series 

We provide two different examples where conditional expectations can be calculated by Dyson series. The 
first example is a very well-known example, but it illustrates nicely the computation of Dyson series in case 
the random variable F (seen as a functional of Brownian motion) is not path-dependent. In the second 
example, the functional F is path-dependent. 
Example 1: 
Let t > T and: 



FlU) = Pxp( ; -) 

1 ' V^T yy 2{t-T)> 



We calculate the Dyson series: 



M{t, U ) = M(T, ^H) + i C D 2 -±= exp(- ^ ^ )ds 1 



< - D 2 . ... 

2j t S1 V^T " v 2(r-T) 

1 , W 2 (t,uj), l /m ,. VL 2 (i,w) 1 . , W 2 (t,w). 

i w2 &") , , /, /tin ^ f i , "M , 

ex P(-^ — ^n) + (* - T )l^ / ^ ex P(-777Z — ^t)- + - 



Vr- 


T 


1 




Vr- 


T 


1 




vV- 


— e 

f 



2(t~T)' y 'dT\^7^T ^ y 2(r-T)' 



2(r-*) 

Observation: We deliberately took r > T so that the functional F would be infinitely Malliavin differen- 
tiable. It remains to be seen whether proper convergence results can be obtained when r J, T. 

Example 2: a path-dependent functional 

Let F{uj) — exp(— J Q T W{s,oS)ds). We introduce informally the process W (singular white noise) in order to 
ease the calculation of the Malliavin derivatives. The functional G is indirectly defined as: 

F(u>) = G(W(u)) = exp(- / f W[u,u)duds)) 

Jo Jo 

= exp(- / (T — u)W(u,w)du)) 
Jo 
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Malliavin derivatives can be computed formally as regular derivatives of G. Thus: 

D 2 S F = (T — s) 2 F 
F(w*(w)) = exp(- / W(s,oj)ds-W(t,u))(T-t)) 







The Dyson series becomes: 

M(t,w) = F(w*(w)) + £ / (T - «i) a F(w*(w))dsi + 



1 ' T 



2 



= exp(-y W(a,u)da-W(t,u)(T-t))(l+^J {T - Sl ) 2 d Sl ) + ..) 

/"* 1 

= exp(- / W / (s,w)ds)exp(-VF(i,w)(T-f) + -(T-t) 3 ) 

Jo 6 

3 Representation of Solutions of Path-dependent Partial Differen- 
tial Equations 

We now introduce some key concepts of the functional Ito calculus introduced by Dupire. For more in- 
formation, the reader is sent to Cont and Fournie (2011), which we copy herefater almost verbatim. Let 
X:[0,T] x fl — > R be a continuous semimartingale. The paths of X lie in Co([0, T], R). We also consider the 
space Z?([0,i],R) of cadlag functions with values in R. For a path x £ £)([0,T],R) we denote by x(t) the 
value of x at i, and by x P — (x(u),0 < u < t) the restriction of x to [0,i]. For a process X, we denote by 
X(t) the value of X at t, and by Xf = (X(u),0 <u<t) its path on [0,t]. Consider a path xf e L>([0,i],R). 
For T > t we define the horizontal extension x P h E Z3([0,T],R) of by: 

p . . _ J x(m) if u e [0,t] 

x t,hW ~ \ x (t) if ue]t,t + h] 

A non-anticipative functional F is a family of functionals F = (F(., £)) tg [ 0] T] where: 

F(.,t) : D([0,i],R) -» R 

x F(a;,*) (8) 

is measurable with respect to Bt, the canonical filtration on D([0, t], R). A non-anticipative functional is said 
to be continuous at a; £ D([0,T],R) if: 

Ve > 0,3r] > 0,Vx' then d^x.x') < n =>• |F(x,t) - F(x',t)\ < e 

For ft 6 R we define the vertical perturbation xf'^of xf as the cadlag path obtained by shifting the endpoint 
by h: 

xf' (u) — x\ (u) + hl[t = u] 
The horizontal derivative at i 6 D([0,T],R) of a non-anticipative functional F = (F(.,t)) is defined a£]: 

F{x p ,t + h)-F(x p ,t) 
A t F(x)= lim ± Llli (9) 

whenever it exists. The vertical derivative of F(.,t) at x is defined by: 

7i->0+ ft 



1 We depart here from the notation in Cont and Fournie to use the original notation in Dupire. The reason is the risk of 
typographical confusion betwen the horizontal derivative and the Malliavin derivative. 
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We define C 1,fe ([0,T]) as the set of functionals which are left-continuous, horizontally differentiable with 
A t F continuous at fixed times, and k times vertically differentiable with A^F left-continuous. We define 
Cj' fe ([0,T]) as the set of functionals F e C 1 ' fc ([0,T]) such that A t F, A X F, A k x F e B([0,T]), the set of 
non-anticipative functionals which satisfy definition 2.4 in Cont and Fournie (2011). 

THEOREM 4 Suppose that a,b and G are non-anticipative C£'°°([0, T]) functionals. Suppose there exists 
a unique solution v <G C^ ,oo ([0, T]) of the semilinear PPDE 

A t v + a(x P ,t)A x v + ^b(x P ,t)A xx v = 

v{x p T ,T) = G{x p T ,T) 

Let X t be a stochastic process satisfying: 

dX(t) = a{X P , t)dt + b{X P , t)dW(t) 
Suppose G(X p (lo),T) = F(u) for some F e Then, in L 2 {P): 

v{t,X p (uj))=TeM\f t D 2 s ds)F{u\u)) 

Proof 

The smoothness conditions on a, b, and G imply the smoothness conditions (HI) and (H2) in Peng and 
Wang (2011). By their functional Feynman-Kac formula (see also Dupire (2009)) we have: 

v(X p ,t) = M(t) 

where 

M{t) = E[G(X p ,T)\F t ] 
By proposition 5.6 in Oksendal (1997), for any I and t G [0, T] we have: 

A l x v(X p ,t) = D l t M(t) 

4 Conclusion and Future Work 

The main result of this paper is theorem 3. Formally it contributes to building yet another (small) bridge 
between stochastic processes and quantum field theory. For future work, we intend to design and analyze 
new numerical schemes that implement the Dyson series to solve BSDEs. The main weakness of theorem 3 
is that it currently requires the functional F to be infinitely Malliavin differentiable. A lot of work seems 
to be needed to relax some of these smoothness requirements. In the Markovian case (i.e., F(w) being a 
function of W(T, to)), it is known that theorem 4 generally applies when F is only twice differentiable, as can 
be proved from the theory of quasilinear parabolic PDEs (see for instance McOwen (1995), chapter 11). For 
path-dependent functionals, the picture is much less clear. 

Theorem 3 can certainly be extended to a filtration generated by several Brownian motions, and probably 
to Levy processes. A generalization from representation of martingales to representation of semimartingales 
would also be interesting. Because of the growing importance of fractional Brownian motion, an extension 
of this theorem to fractional Brownian motion would be desirable for applications. 
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5 Appendix 

Following Oksendal (1997), we say that a real function g : [0,T] — > R™ is symmetric if: 
for all permutations a of (1, 2, .., n). If in addition 

ll„l|2 / 2 



lbllL2([o,T])» = / 5 (xi,..,x n )dxi..dx n < oo 

[0,T]« 

Then we say that g e L 2 ([0,T] Tl ) the space of symmetric square integrable functions on [0,T] n . The Wiener 
chaos expansion of F: 

oo 
m=0 

where {f m }m=o i s a unique sequence of deterministic functions in L 2 ([0,T] n ), and 

im(/m) = / [ m •■ [' f(t 1 ,..,t m )dW(t 1 )dW(t 2 )..dW(t m )i{m>0 



Jo Jo Jo 
Wo) = fo 

T 

A notation we use for the Skorohod integral of H is J H(s)SW(s). 

o 

5.1 Proof of Theorem 1 

Let t — mA and T = (m + 1)A. We remind the reader once again of proposition 5.6 in Oksendal(1997), 
namely that, if F £ Oi 2 and ^[FlJ^] <E ©1,2 (see Oksendal fro a definition Di 2) then, for i < s 

D t {E[F\Fs\) = E[D t F\F s ] (10) 

Let t + — t + e for a small e. Using (fT0| and the Clark-Ocone formula (see, e.g. Nualart (1995)), we get, for 
t + < T and any u: 



r t+ 

E[D[F\F t+ ] = E[D l u F]+ E[D s E[D l u F\F t+ ]\F s ]dW(s) 

Jo 

= E[D l u F]+ E[E[D s D l u F\F t+ ]\F s }dW(s) 
Jo 



(11) 
(12) 

= E[D l u F] + f E[D s D l u F\F s ]dW(s) (13) 
Jo 

= E[D l u F] + [ E[D s D l u F\F s ]dW(s) (14) 
Jo 

- [ E[D s D l u F\F s ]dW{s) (15) 
Jt+ 

= E[D l u F\F T ] - [ E[D s D l u F\F s ]dW(s) (16) 
Jt+ 

Since F is assumed discrete, for s € (t, T] we have: 

D s D l u F = D T D l u F (17) 
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Thus 



E[D l u F\F t+ ] = E[D l u F\F T ] - J E[D T D l u F\F s ]dW(s) (18) 

Taking u = T, we obtain: 



E[F\F t+ ] = E[F\Tt) - I E[D l T F\T Sl ]dW{ Sl ) (19) 

Jt+ 

= E[F\F T ] - [ f E[D X T F\F T ] - f E[D 2 T F\F S2 ]dW( S2 )SW{ Sl )] (20) 

Jt+ Js, 



where (|19|) follows from (|18JI with I — 0, and (|20[) follows from (|18|) with / = 1. We continue the expansion 
iteratively until we calculate the L th Malliavin derivative, after which all terms are zero. Then, by continuity 
of martingales generated by Brownian motion, we conclude that: 

E[F\T t ] = E[F\Tt\- E[D l T F\T T ]+ E[D 2 T F\F T ] - (21) 
.. + {-l) L [ E[D%F\F v ]5W{s L )..8W{s x ) (22) 



where these integrals are iterated Skorohod integrals. For convenience, we define the iterated time/Skorohod 
integral N(b, so). Let b be a binary vector of dimension r (which will be clear from context). We define —b 
as the same vector without the first component. It is thus a vector of dimension r — 1. For instance, if r = 5 
and: 

b=[l 1 ]=>•-&= [0 1 0] (23) 
The iterated time/Skorohod integral is: 

f 1 r = 

N(b, t) = \ jf N(-b, s)ds if r > 1 6 r = 1 (24) 
I jf JV(— b, s)<W(s) r > 1 b r = 

For instance, with 6 as in (|2"3"|) we have: 

W(M) = J J J J I SW{s 5 )5W(s A )ds 3 SW(s 2 )ds 1 

t Si S2 S3 S4 

We define 

Af,, (t) = S T [£>^F|J- T ] 

T T T 

Mj, 9 (t) = J .. J J E[D l T F\F T }SW( Sg )SW( Sg ^)..SW( Sl ) g>0 

t S„_2 S„_i 



Lemma 1.1 



S 



min(/,L — I) 

M u (t)= E E[D l + h F\F T ] J2 

h=0 b 1 + ..+b l =h 

Proof of lemma 1.1 

From the formula for the Skorohod of a process multiplied by a random variable ((1.49) in Nualart (1995)), 
we calculate: 



Mi,i(«,_i) 



E[D l T F\F T ]5W{si) 



E[D l T F\F T }N(0, aj.O - E^ 1 F\F T ]N {I , s^) if I < L 
E[D l T F\F T ]N(0, sj_i) if I > L 



We suppose by induction that: 



MiJsi- 9 )= £ E[D l + h F\F T ] Y, (-l) k N(b,s^ g ) 

h=0 bi+..+b„ = h 



Reapplying (1.49) in Nualart (1995), we obtain: 



Ml,g +1 (Sl- g -l) 



min(p,n— I) 



<T E[D l + h F\F T ] f 
E[D l + h+1 F\F T ] / 

J s l-g-l 



Y (-l) h N(b, Sl - g ) 
b l +..+b g =h 

J2 (-l) h N(b, Sl . g ) 



b 1 + ..+b„=h 



SW( Sl . g ) 



dsi- 



H-..-|-0<,= 

min((/+l,n — I) 

]T E[Dp h F\F T ] (-l) h N(b,si- 9 -i) 

h=0 b 1 + ..+b g+1 = h 



(25) 



With lemma 1.1, (|21j) results in: 



E[F\F t ] = £(-l)<M M (i) 

1=0 

L min(Z,L — V) 

= E E ^tI^^t] e (-^(m) 

1=0 h=0 b 1 + ..+b,=h 

L 

= Y,l{™,l)E[D l T F\F T ] 

1=0 



(26) 
(27) 
(28) 



where 7(771, Z) does not depend on i 7, or L. One possibility to calculate 7(771, Z) is to use lemma 1.1. For 
I < L/2 we have for instance: 



7(771, Z) = ( 



!)E E (-^(m) 



/i=0 bi + ..+b,=/i 
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This is rather complicated. However, since (|28|) holds for any diffcrcntiablc F, we resort to a simpler strategy. 
Our strategy is to vary F to determine recursively 7(771, 1). For simplicity, we take m = 0, i.e., t = 0. Clearly 
the first coefficient (take F=constant) is: 

7(0,0) = 1 

To determine the second coefficient, 7(0, L), with L = 1, we choose a function F such that D L+1 F = 0. The 
only such function is a linear function of W(T). We thus put in (f28|) F = W(T) and calculate: 

1 

E o [F]^-y{0,l)D l T F 
1=0 

In other terms: 

= 7(0, 0)W(T) + 7(0, 1) * 1 

Thus 

7 (0,L) = -W(T) 
The general structure of the recursion is then: 

s E\F] - Y^" 1 7(0, l)D l T F , N 

7(0, L) = 1 J pip T (29) 

Clearly, formula (|29|) applies for any coefficient L. To calculate 7(0, L) for L = 2 we thus pick F = W 2 (T) 
(so that D L+1 F = 0) and obtain: 

e{w 2 (t)}-EL j(o,i)d 1 t w 2 (t) 



7(0,2) 



£> 2 W 2 (T) 
T - W 2 (T) + W(T) * 2W(T) 
2 

W 2 (T) +T 



2 

We complete the proof by induction. Suppose that, for 



7(0, L-1) = 1{L even} - E7 (Q,,W±±_^ (30) 



(L/2)! £j '^'(L-i-l)! 



We select F = W L (T), which we insert together with §5U§ in ([25) to arrive at 

7(0, L) = 1{L even} - ^^.Ot^T < 31 > 



i=0 



The formula for general m obtains by replacing H^(T) by W((m + 1)A) — W(mA), that is choosing for test 
functions F above the successive powers of W((m + 1)A) — W(mA). 
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5.2 Proof of Theorem 2 



Let t — mA and T = (rn + 1)A. By theorem 1: 

T T 

E[F\F t ] = F% + (~1) L j .. J E[D^ 1 F\T SL+1 ]6W(s L+1 )...6W(s 1 ) (32) 

si=t S L + 1 =S L 

We want to reverse the order of integration in (|32[) . and thus show an elementary Fubini-type theorem for 
Skorohod integrals from first principles. Let /„ be a deterministic function of n + 2 variables. Suppose the 
function is symmetric in the first n arguments. We define the following symmctrization operators: 



Sym n ' pen (f n ){t u ..,t n ,v,u) = -^—[f n (t 1 ,...,t n ,v,u) + 

n + 1 



/n(*l) V,t n , u) + .. + f n (v,t 2 ..,t n ,ti,u)} 



Sym n > last (f n )(t u ..,t n ,v,u) = — — -[f n (t 1 ,...,t n ,v,u) + 



n + 1 

/„(*!, ...,t„_i, U, U,t n ) + .. + f n (u,t 2 , ..,t„,V,ti)] 



Sym n+1 *™{f n )(t u ..,t n ,v,u) = -^—:[fn{ti,...,t n ,t n+l ,v) + 



n + 2' 

fn{t\, f„, V, f n +i) + .. + /n(f, t2, ■•, ^n, 1 1 , tn+l )] 

5ym n+1 ' , °'*(/„)(ti ) ..,t Tl ,«,«) = — — - [f n (t u ...,t n ,t n+1 ,u) + 

n + 2 

f n (tl, i n , It, + ■• + /n(u, tl) •■, t n , t n+ i,ti)} 

The operator Sym n ' pen (f) sums all the permutations of any of first n variables with the penultimate vari- 
able ("u"). The operator Sym n,last (f) sums all the permutations of any of first n variables with the last 
variable("u"). The operator Sym n+1 ' pen (f) sums all the permutations of the penultimate variable ("v") with 
any of the other n+1 variables. The operator Sym n+1:last (f) sums all the permutations of the last variable 
("it") with any of the other n + 1 variables. The following elementary lemma translates into our notation the 
fact that the operations of symmetrization with respect to the penultimate variable (" v" ) and symmetrization 
with respect to the last variable ("u") are commutative. So, without proof, we state: 
Lemma 2.1 

Sym n+1 ' pen (Sym n ' last (f)) = Sym n+1 > last (Sym n ' pen (fj) 
Of course, lemma 2.1 applies to 

g(ti, t ni v, u) = t n , v,u)*l[v> u] 

Lemma 2.2 

Suppose H(s) is an J^-adapted process. Then: 

H(s 2 )dW(s 2 )6W{si) = j { H(s 2 )dW(s 1 )dW(s 1 ) (33) 

si—O J S2—S1 J S2—O ^ si—0 

Proof of lemma 2.2 

Since H(s) is adapted, the iterated Ito integrals on the right of (|33p are identical to Skorohod integrals, and 
similarly for the inner integral on the left. It is then sufficient to prove: 

j H{s 2 )8W{s 2 )5W{s x ) = f f 2 H(s 2 )SW(s 1 )SW(s 2 ) (34) 

Si=0 J S2—S1 J S2—O J Si— 

Let: 
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f H{s 2 )8W{s 2 ) = f; J n+ i(^(, s 2 , si)) 

•I S 2 =8i n 



'S2=S1 „ = 

where, for some /„ symmetric in the first n variables: 



and 



where 



h l n (t l ,...,t n ,s 2l s 1 ) = Sym n < pen (f n (t 1 ,...,t n ,s 2 ,s 1 ) * l[s 2 > si]) 

pT pT oo 

/ / H( S2 )SW(s 2 )5W( Sl ) = W»»(-> s x )) 



We see that: 



where 



Then 



where 



ft"(ti, t B , s 2 , si) = Sym n+1 ' lo '*(/i 1 (ti, t n , s 2 , ai)) 
ff(« 2 )W(*l) = X)l„+i (<£(., «2,*l)) 

-1=0 n =0 

g*(t 1 ,...,t n ,s 2 ,s 1 ) = Sym n ' last (f n (t 1 ,...,t n ,s 2 ,s 1 ) * l[s 2 > s{\) 

pT ps 2 00 

/ / H( S2 )6W(si)6W( Sl ) = J2ln +2 (gU,S2,-)) 

J s 2 =0 J s 1= „_q 



g 2 n (h, t n , s 3 , Sl ) = 52/m" +1 ' pe "(^(ii, -, *», *a, 
The left handside of f[M|) is then equal to: 

oo 

I n +2(Sym n+1 ' last (Sym n ^ en (f n (h, t n , s 2 , fll ) * l[s 2 > Sl ]))) 
The right handside of is equal to: 

oo 

I n +2(Svm n+1 ' pen (Sym n ' last (f n (h, t n , s 2 , Sl ) * l[s 2 > Sl ]))) 
By lemma 2.1, (|35|) is thus equal to (|36|) . ■ 



Suppose if (s) is an ^-adapted process. Applying lemma 2.3 we obtain: 

/ .. f H(s n+1 )SW(s n+1 )..6W( Sl ) = [ .. [ 2 H(s 2 )5W( Sl )..SW(s 

J S±— t J Sn-^i—Sn J S n _|_l— t J S\—0 

With the help of ([37]), with E[D^ +1 F\F S ] replacing H(s), relationship fl32J} becomes: 

Ft] = F r L n + {-l) L J .. J E[D^+ 1 F\F SL ]dW(s 1 )...dW( SL+1 ) 

Sl=t S L + 1 =t 

By lemma 5.7.2 in Kloeden and Platen (1992), we have: 
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E[{E[F\F t ]-F^f\F t ]< sup E[{E[D^ F\F s ]) 2 \F t ]——- (39) 

s£(t,T\ { L + 1 ) 1 

For s G (t,T], the following process is a martingale: 

X( S ) := £[y^ +1 ^| F s ] (40) 

Thus: 



sup ^[X( S ) 2 |J" t ] = sup E[X{t) 2 + f d[X,X](u)\T t ] 

se[t,T] se[t,T] Jt 

E[X(t) 2 + J d[X,X]{u)\F t ] 
E[X{tf + j\E u {D L u + 2 F]fdu\T t ] 



< E[X(t) 2 \T t ] + A sup E[D L U +2 F\F U ] 2 
ue[t.T] 



By induction, 



sup E[X(s) 2 \F t ] = J2^(E[Dt +l F\T t ]) 2 

»e[t,T] l=0 

By assumption, all Malliavin derivatives are bounded. Thus, if A < 1, there is a constant M so that, for all 
I: 

sup E[X{s) 2 \F t ] < -*L- (41) 

«€[t,T] 1 - A 

By putting together (J3SJ, (gDJ), and ([ST]), we have: 

Ai+l TUT 

E[{E[F\F t ] - FZf\F t ] < (L + 1)!1 _ A ( 42 ) 

5.3 Proof of Theorem 4 

Let _F( T,n ) be the approximation "by simple functions" of the random variable F. More precisely, we define 
the discretizing map t( T '™) 



t if t<T-l/n 

T-l/n if T -l/n<t <T 



For each / OT we define the approximation by simple functions fm 

f£> n Hh,t 2 ,..,t m ) = f m {r {T ' n \tx),T^ n \t 2 ),..M T ' n) {t m )) m>l 

AT,n) _ , 
Jo — JO 

We now define: 

oo 

F™ = £ I m (/£») 

m=0 

We also define J-Q{T, n ) to be the sigma-algebra generated by Brownian motion under the condition that: 
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W(s) = W(T- l/n) T -l/n<s <T 

Lemma 4.1 

E[(F( T -")-F) 2 |jr G(T , n) ]=0 
Proof By orthogonality of the Wiener chaos expansion, we need to prove that 

E[{Im{f^ n) -fm)f\r G{T ^]=Q 

The case m = 1 spells: 

(rp \ ( T 1/iV ^ 

^ A (T ' n) (*i)-/i(tiW*i)J l^ 0(TlB) ] = e[U ' /£ T - n HM -/*>(*! 

= o 

Let 

x^(t) = f f m 1 •• f 2 f£' n Xh,..,t m - 1 ,t)dW(t 1 )dW(t 2 )..dW(t m - 1 ) 

Jo Jo Jo 

X m -i{t) = [ f m 1 .. f' f m (t 1 ,..,t m - 1 ,t)dW(t 1 )dW(t 2 )..dW(t m - 1 ) 
Jo Jo Jo 



2 



By induction: 



Thus 



X ( ^(t)=X m _ 1 (t) for t<T- l/n 

E[{I m {fS' n) )-Im{fm)f\T G ^ n) ] = 
E{(J o X^(t m ) - X m ^(t m )dW(t m )) 2 \T G(T ^} = 

E[( X^(t m )-X m ^{t m )dW{t m )) 2 ] = 

Jo 



Lemma 4.2 For any I integer: 

E[(D l t F™ - D l t F) 2 \T GlT , n) } = 

Proof 

We only prove the case I = 1. The case of higher derivatives is similar. By lemma 4.16 in Oksendal (1997): 

oo 

D t F= J2 mI m-l(fm(;t)) 
m—1 

and also 



m—1 
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Then 



E[(I m -!{f^%,t) - fm(.,t))f\F G ™} = 

E[( I X™(t m ) X m ^{t m )dW(t m )f\F G{T ^] = 
Jo 

fT-1/n 

E[( / X™{t m ) - X m ^{t m )dW{t m )) 2 ] = 
Jo 

The convergence of D t F^ T ^ to D t F follows by orthogonality of the Wiener chaos expansion. 



Lemma 4.3 
Proof 

By proposition 5.5 in Okscndal (1997) 



E[F (T,n) _ F \JT T _ 1/n ] = 



where: 



By construction: 



Thus: 



Lemma 4.4 



E[Im(fm)\^T-l/n] — Im{f m lfF_i/ n j 



C/ml^"])(il7 ■;t m ) — / m (<l, ••,im)l[0,T-l/n](^l)---l[0,T-l/n]( i m) 



(fm ' ^[0,t] /ml^ jt ])(il) ■•) im) — 
^ m (/i T ' n) -/ m )I^T-l/„]=0 



lim £;[(n( J B[^ T '")|J- T _ 1/n ]-^ T '"))-^|F) 2 |J- G(T , ri) ]=0 



Proof 

By the backward Taylor expansion (theorem 3) for any n: 



E[F™ |JT T _ 1/7i] = £ 7(T _ 1/n) l )D l TF (T,n) 



1=0 



We remind the reader that along the path lj t 1 /™(cj) we have: 

W{s,uj T - 1 ' n {u)) = W(T-l/n,u) T — 1/n < s <T 

Thus: 



7 (T- 


1, w T 


" 1/n (w)) 


7 (T- 


l/n,2,w T ~ 


- 1/n (u)) 


7 (T- 


l/n,3,w T ~ 


- 1/n {u:)) 


7 (T- 


1/n, 4, 


- 1/n M) 


7 (T- 


1/n, Z,w T ~ 







1 

2^ 


1 

0(-U for / > 4 

n 6 
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By the Cauchy-Schwartz inequality, for I > 4: 



= °(4) 

Thus: 

E[E[F™\F T _ 1/n ] 2 ] = E[(F™ + ^D l T F^ + 0{ \)f\F G ^] 
and the conclusion follows. ■ 

Lemma 4.5 

lim E[{n{E[F\T T -i/ n ] — F) — ^F) 2 |F g( t,„>] = 

Proof 

We decompose the integrand into: 



-V 

with: 



(F[F|F T _ 1/n ] — F) — \d 2 t F = A n + n(F - F< T -">) 
A n - n(F[F|F T _ 1/n ] - F™) - ±D 2 t F 



2 

By the Cauchy-Schwartz inequality and lemma 4.1, 



lim E[nA n (F - F^)\F G(T , n) ] < lim ^/ E[Al\F G{T , n) ]J E[v?{F - F^)y\F &T , n) ] 



n— >oo 



lim y /E[A n \F a <.T,n->] Jim v/F[n 2 (F - F^n)^ 




n— »-oo n— >oo v 



Applying lemma 4.1 and developing (|45p , we have: 

hm F[(n(F[F|F T -i/„] - F) - ^F) 2 |Fg(^»>] = hm F[^ 2 J F G( r,„)] 
We further decompose A n into: 

A„ = B„ + i(£ 2 F( T < n )-i} 2 F) 

with 

B n = n{E[F\F T _ l/n ] - F<*») - ^D 2 F (T ' n) 
Again, by the Cauchy-Schwartz inequality and lemma 4.2: 

lim E[B n (D 2 F^ - D 2 t F)\F G{T , n) ] = 



Applying lemma 4.2 and developing (|4"T)) . we have: 

lim FL4 2 |F G(T ,„>] = lim F[B 2 |F G(T ,„>] 
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We further decompose B n into: 



B n =C n + nE[F- \F T -i/n] (49) 



where: 



Thus, by lemma 4.3: 



C n = n(E[F<?>»\F T _ 1/n ] - F™) - \d^F^ 
hm C n E[F-F™\F T ^ 1/n ]=0 



Thus, by lemma 4.4: 

lim E[B 2 n \F G{T , n) ] = lim E[C 2 \F G(T , n) ] (50) 

lim E[(n(E[F™\F T . 1/n ] F™) - l -D 2 F^) 2 \F G{T , n) ] 

o 



1/nj - ' 



Putting together (|46l) . (f4"5)) and (|5"0")) . the lemma obtains. ■ 

We now replace T by s and F by E[F\F S ] in lemma 4.5 to conclude that, for any s G [t,T]: 



hm E[{n{E[F\F^ 1/n ] - E[F\F S }) - T;D 2 E[F\F s }) 2 \F G(s . n) } = (51) 



1 

2 

We define F H ( S ,t) to be the sigma-algebra generated by (for all w G f2). We clearly have 

1 



lim E[(n(E[F\F s _ 1/n ] - E[F\F S }) - -D 2 s E[F\F s }) 2 \F H{s - 1/n . T) ] = (52) 



We write: 



M(a,u) = E[F\F.](u) (53) 
M(s,ui t (ui)) is J^-measurable for each s G [t, T]. Thus, we apply Fatou's lemma to (|52|) . so that, for any 

lim r» 2 (M(»,w'- 1 / n (w)) - M(a,u'-^ n (u))) - )-D 2 s M{s,u: s ~ 1/n {uj)) =0 in L 2 {P) 
n— >oo 2 

From now on, all equality statements for random variables are understood to be in L 2 (P). A fortiori, the 
sigma-algebra generated by w'(u;) (for all oj G O) is included in F H ( S -i/ n ,T) , for any s G [£+— ,T]. Thus, with 
t<s~At: 

M(s-AWM)-M(s^)) = i^ M{Wiu)) (54) 
At^o Ai 2 v / / v y 

M(T,w*(w)) = F(u*(u;)) (55) 

The remainder of the proof is standard, but we provide it nevertheless for completeness. Let T-Lt be the Banach 
space of linear operators mapping a bounded Ft -measurable random variable to a bounded Ft- measurable 
random variable in L 2 (P), with the operator norm. We now define the propagator P(s, r) G Ti s as a linear 
operator mapping M(r,w t (w)) to M(s,w*(a;)) in £ 2 (P) with A/ defined by ([55]). with t < s < t < T. We 
have: 

IIPftrJik-^MW*-"! 



£[M 2 (£)|.F ff(t .T)] 

Thus F(t, t) is bounded. Also, the semi-group property obtains. For t < s < t < T: 
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P(s,t)P(t,T) = P(s,T) 

Equation (|5^|) and (|55|) become then: 



(56) 



d -^l = \dIP{s,T) t<s<T (57) 
P(T,T) = I (58) 



The evolution equation (|57|) and (|58|) is equivalent to the integral equation: 

P(t,T) = I+±J D 2 P(s,T)ds 
We define the linear operator F(£, r) : — > % by: 

VQt,T)P(t,T) = I+±J D 2 s P(s,T)ds 
We now show that the fixed point equation: 

P(t,T) = V(t,T)P(t,T) (59) 

has a unique solution. We show that V^(T — (5, T) is a contraction on the space Ht-8 for sufficiently small 5. 
We define the operator norm 

V(t,T)| = SUp 



We now estimate: 



II^T-^TJti-yCT-tf.TjwHoo : = sup \\V(s,T)u(s,T) — V(s,T)v(s,T)\\ 

T-8<s<T 

< isup f \\D 2 s u(s,T) - D 2 v(s,T)\\ Hs ds 

1 Jt-8 

By assumption the Malliavin derivatives are bounded, i.e. there is a constant C such that: 

\\D 2 s u(s,T) - D 2 v(s,T)\\ n3 < C\\u(s,T) - v(a,T)\\ Hm 
We conclude that (with obvious definition of the norm on the right handside): 

\\V(T - S,T)u - V(T - «,T)t;||oo < CS\\u - u||oo 

It follows that, if 6 < 1/C then V(T— 5, T) is a contraction on the space of operators Ht-8 ■ By the contraction 
mapping theorem there is a unique solution to the equation P(T — S,T) — V(T — 5, T)P(T — 6, T) and the 
solution can be constructed by the Picard iterations: 



P (T-S,T) = I (60) 
P n+1 (T-S,T) = V(T-8)P n (T-S,T) 



Following verbatim the argument in Zeidler (2006) p. 388 and 389: 
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P(T-5,T) 

By the semi-group property (I56P : 

P(T-2S,T) 

By covering the interval [t, T] by overlapping intervals of length less than 1/C, theorem 4 obtains. 
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/ + E/ -/ (h k T(D 2 Ti ,..Dl k )dr,..dr k 
k=1 Jt-s Jt~s a 

T(exp(i f T D l ds) 

L JT-8 



= P(T-2S,T-5)P(T-5,T) 

= T(exp(i f 5 Dlds)T(eM\ f D 2 s ds) 

L JT-2& 1 JT-S 

= T(exp(i f Dtds) 

L JT-25 
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